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Abstract 



We study supersymmetric extensions of the Standard Model with small R-parity and 
lepton number violating couplings which are naturally consistent with primordial 
b^i nucleosynthesis, thermal leptogenesis and gravitino dark matter. We consider super- 

,^ I gravity models where the gravitino is the lightest superparticle followed by a bino-like 

next-to-lightest superparticle (NLSP). Extending previous work we investigate in de- 
tail the sensitivity of LHC experiments to the R-parity breaking parameter (" for 
^ , various gluino and squark masses. We perform a simulation of signal and back- 

^vq I ground events for the generic detector DELPHES for which we implement the finite 

On ' NLSP decay length. We find that for gluino and squark masses accessible at the 

^-^ . LHC, values of ^ can be probed which are one to two orders of magnitude smaller 

^ I than the present upper bound obtained from astrophysics and cosmology. 



X 



1 Introduction 

Supersymmetric extensions of the Standard Model with broken R-parity have a rich phe- 
nomenology [1-5]. In most models rather large R-parity violating couplings are considered, 
which lead to prompt decays of the lightest superparticle in the detector. In models where 
small R-parity violating interactions generate neutrino masses, macroscopic decay lengths 
up to 1mm are obtained [6]. In the case of gauge mediated supersymmetry breaking, 
R-parity violating decays then compete with R-parity conserving decays where the final 
state contains a gravitino [7]. 

If the gravitino is the lightest superparticle, its decays are doubly suppressed by the 
Planck mass and the small R-parity breaking parameter. Hence, its lifetime typically ex- 
ceeds the age of the universe by many orders of magnitude, so that it remains a viable dark 
matter candidate [8] . In the case of very small R-parity breaking couplings, as they occur 
if R-parity is spontaneously broken at the grand unification scale, one obtains a consistent 
cosmology including primordial nucleosynthesis, thermal leptogenesis and gravitino dark 
matter [9]. In the following we shall study the collider phenomenology of these models 
extending our previous work [10]. 

Decaying gravitino dark matter leads to a diffuse gamma-ray flux [8,9, 11]. For a bino- 
like NLSP the matrix elements for gravitino decay and NLSP decay are directly related [10]. 
Together with the lower bound on the gravitino lifetime, and the corresponding upper 
bound on the R-parity breaking parameter (, which is derived from the diffuse gamma-ray 
flux observed by the Fermi-LAT collaboration [12,13], one obtains a lower bound on the 
NLSP decay length. The estimates in [10] led to C ^ 3 x 10~^, and in the more detailed 
analysis of gamma-ray lines in [14] a slightly stronger bound was found, C ^ 2 x 10~^. 
For gravitino masses 0(100 GeV) the upper bound becomes more stringent: ( ~ 10^^ [10]. 
This range in the parameter ( corresponds to NLSP decay lengths varying from (9(50 cm) 
to 0(500 m). 

Large macroscopic decay lengths are of great help in the search for decaying NLSPs. 
This remains true if the decay length is larger than the size of the detector since a sizeable 
fraction of NLSPs may still decay inside the detector. This has been studied for neutral 
[15] as well as charged [16] NLSPs. Neutralino decay lengths varying from 0.1mm to 
100 m also arise in models with generalized gauge mediated supersymmetry breaking [17]. 
Alternatively, charged [18] and neutral [19] NLSP decays have been studied for models 
where the decay lengths are so small that no displaced vertices are observed and R-parity 
breaking Yukawa couplings determine the hierarchy of decay channels. In this case multi- 
lepton events, and their flavour structure, are of crucial importance. 

The subject of this paper is a quantitative analysis of neutralino NLSP decays at the 
LHC in the case of very small R-parity breaking, ( < 3x 10~^. The goal is the determination 
of the sensitivity in ( for varying gluino and squark masses. To achieve this we perform a 
simulation for the generic detector DELPHES [20] of signal and background events. We focus 
on events with a clean signature: cascade processes with jets where one of the produced 
neutralino NLSPs decays into Z-boson and neutrino, with the subsequent decay of the 
Z-boson into a muon pair. This allows us to determine a conservative 5a discovery range. 



Finally, we estimate the discovery reach of the LHC if all NLSP decays are taken into 
account. 

The paper is organized as follows. In Section 2 we recall the main ingredients of the 
considered model of bilinear R-parity breaking, as well as the dominant NLSP decays 
into gauge bosons and leptons. In Section 3 we discuss qualitative signatures of the events: 
NLSP (3'~f distribution, p spectrum, and the number of final state leptons. In Section 4, five 
representative benchmark points are defined and the simulation of signal and background 
events is described with emphasis on the muon reconstruction. The numerical results for 
the chosen NLSP decay channel Z{n'^n~)i> are given in Section 5, together with estimates 
for the discovery reach if all NLSP decays are included. We conclude with a summary in 
Section 6. 



2 Connecting neutralino and gravitino decays 

2.1 Bilinear R-parity breaking 

We consider supersymmetric extensions of the Standard Model with bilinear R-parity 
breaking (cf. [4,5]) as they are obtained if the spontaneous breaking of B— L, the difference 
of baryon and lepton number, is related to the spontaneous breaking of R-parity [9] . Mass 
mixing terms between lepton and Higgs fields then appear in the superpotential^, 

AW = fi,HJ, , (1) 

as well as the scalar potential induced by supersymmetry breaking, 

- A£ = B.HJi + mU\Hd + h.c. . (2) 

These mixing terms, together with the R-parity conserving superpotential 

W = fiH^H, + h^l^q^u'^.H^ + hf^d'rq.Ha + hy.e'jH, , (3) 

the scalar mass terms 

-Cm = mlHlH^ + m^i/ti/, + {BH^H, + h.c.) 

+ ^IM + mlfife'l + mjgjgi + mliufu^i + ml^fdl , (4) 

and the standard SU(3) x SU(2) x U(1)y gauge interactions define the supersymmetric 
standard model with bilinear R-parity breaking. Note that the Higgs mass terms m\ and 
m^ contain the contributions both from the superpotential (3) and the soft supersymmetry 
breaking terms. For simplicity, we have assumed fiavour diagonal mass matrices in (4). 



^Our notation for Higgs and matter superfields, scalars and left-handed fcrmions reads: Hu = {Hu, hu), 
H — \H 1 H ) etc. 



As discussed in [10], it is convenient to work in a basis of SU(2) doublets where the mass 
mixings /ij, Bi and mf^ in Eqs. (1) and (2) are traded for R-parity breaking Yukawa cou- 
phngs. This can be achieved by field redefinitions: the standard rotation of the superfields 
Hd and k, 

Hd = H'^- ^il'i , h = l'i + tiH'^ , e, = — , (5) 

/i 

followed by a non-supersymmetric rotation involving all scalar SU(2) doublets, 

H'd = H'd - e-/-' , eHl = eH'* - e'-l" , /• = /■' + e[H'^ + e'-eH'^ , (6) 

where e^ and e'- are functions of B, B^, mj^, mf^, m^ and m\ [10]. 

The R-parity breaking Yukawa terms contain couplings between gauginos, lepton dou- 
blets and Higgs doublets. After electroweak symmetry breaking, {H^ = f„, {H^) = Vd, 
one obtains new mass mixings between higgsinos, gauginos and leptons^, 

-ACm D m%—e']hd - mzs^Quib + mzCu,C*ViW^ + h.c. , (7) 

where we have defined 

V Vd Cfj 

rriij = h.jVd , mz = ^ , s^ = ^^, ^ ^,, = ^1 - 4 . (9) 

Here g, g' and h^j are the SU(2) and the U(1)y gauge couplings and the charged lepton 
Yukawa couplings, respectively. 

The diagonal mass terms together with the mixing terms in Eq. (7) represent the 7x7 
neutralino mass matrix of the gauginos b, w^, the higgsinos h^, h^ and the three neutrinos 
z/j, and also the 5x5 chargino mass matrix of gaugino, higgsino and charged leptons. Both 
mass matrices have to be diagonalized to obtain the CKM-type matrix elements of neutral 
{V^^'"^), charged {V^^''^^) and supercurrents {U^^'"''). For the standard supergravity mass 
spectrum with a bino-like NLSP Xi one obtains the R-parity breaking matrix elements 
(s2/3 = 2s^c^) [10]: 

,,r'^-C,i|^(i.o(..f)), ,.0, 

vir>^^(.^i+o(.s,,4)], (11) 



Ml V V A"' 

j,„,,„ _ ^_ ».,..c (Afe-A/i) L^^l ^^m^ ., ,j2, 



mzSy^Cy, (M2 - Ml) / / ml 



Our notation for gauge fields and left-handed gauginos reads: B^, b etc. 



"^ 
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Figure 1: Gravitino decay into photon and neutrino. 



where the photino matrix element is defined as 



(13) 



Note that the charged and neutral current matrix elements agree up to the isospin factor 
at leading order in ml/fi'^, i.e., V^^'-^l = vl^'ll/2. 

2.2 Gravitino and NLSP decays 

The partial width for gravitino decay into photon and neutrino (cf. Fig. 1) is given by [8] 



^S/2M 



3271 



E 



u. 



(7:'^) 



M2 



(14) 



Inserting the matrix element (12) one obtains for the gravitino lifetime to leading order in 

mz/fi [10]: 

-2 



r3/2(7z/) = 1 X 10^^ s 



C 



10 



-7 






100 GeV/ VlOGeV 



m3/2 



where 



c^ = Ecr 



(15) 



(16) 



Based on the Fermi-LAT search for dark matter decaying into two photons, upper bounds 
on the R-parity breaking parameter ( were derived: C ^ 3 ^ 10^^ [10] and C ^ 2 x 10~^ [14]. 
The parameter ( also controls the lifetime of the bino-like NLSP Xi- Fo^' ^ NLSP mass 
larger than 100 GeV decays into charged lepton and VT-boson or neutrino and Z-boson 
(cf. Fig. 2) are dominant [21]. The partial decay widths read 



r{xUZu)=^^Y.K'llffz{m,.)(l + o(s,,^]] ■ (17b) 



2v^' 



fzim^o) ll + O ^s,p-^ 




w^ 
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Figure 2: Neutralino decays into charged lepton and W -boson, and neutrino and Z -boson. 

here V\j^lo ^^^ ^ulo ^^^ ^^^ charged and neutral current matrix elements at leading 
order, which are given in Eqs. (10) and (11), respectively. The function fw,z is a phase 
space factor which becomes important for neutralino masses close to the lower bound of 
100 GeV (cf. Fig. 3), 



fw,zim^o] 



m 



w,z 



Xi 



1 + 2- 



m 



w,z 



X? 



,v 



For large NLSP masses, m^o ^ niz, one has 

BR (x? -^ WH^) c^2 BR(x\^ Z 
whereas in the region m^o c^ 100 GeV 

BR (x? -^ W^l^) -b BR(xI^ Z 
The total neutralino NLSP decay width is given by the sum 

which corresponds to the decay length [10] 



;z/ 



^v) 



Xi 



'^-'^-^■'"^y 100 GeV J 13x10- 



C 



X {2fw{m^o) + fz{m^o 



-1 



-(-f 



(19) 



(20) 



(21) 



(22) 



Note that both, the gravitino and the neutralino NLSP lifetimes are functions just of ( and 
the masses, without any further parameters. This is the case for the standard supergravity 
mass spectrum with a bino-like NLSP. This direct connection between the gravitino and 
NLSP lifetimes is the basis of our analysis. 
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300 



Figure 3: Phase space suppression factors for neutralino decays into W -bosons (red) and Z -bosons 
(black) and leptons. 

3 Decaying neutralino NLSP at the LHC 

In this section we classify the main LHC signatures of decaying neutrahnos. We show that 
for small values of the R-parity breaking parameter C, usual SUSY searches are insufficient 
to find the signal. However, R-parity violation leads to new signals including striking 
secondary vertices at large distances from the primary interaction point. 



3.1 Decay signatures 

Consider for simplicity, the following cascade process: 



19^(19 ^ JJJXiXi , 



(23) 



where g is a squark, ^ is a gluino, and j denotes a jet. The final state neutrahnos decay in a 
secondary vertex into VT-bosons and leptons as well as into Z-bosons and neutrinos. Fig. 4 
shows an example of a decay cascade with muons in final state. The distance between the 
collision point and the secondary vertex depends on the decay width of the neutralino (17) 
and hence on the R-parity breaking parameter (. 

Table 1 summarizes the LHC signatures for sufficiently large values of ( such that it 
is probable that both neutrahnos decay inside of the tracker volume. All the signatures 
contain at least three jets from the antecedent supersymmetric decays. The signatures are 
classified according to the final states in the neutralino decays: leptonic signatures involving 
only leptons in the final state, semi-leptonic signatures involving at least two charged 
leptons and jets, single lepton signatures containing only one lepton, all-hadronic signatures 
where only jets accompanied by neutrinos are present, and finally invisible channels where 
both neutrahnos decay solely to neutrinos. Additionally we single out channels having a 
considerable amount of missing transverse energy ^j- from Z-boson decays, since ^j- is 




Figure 4: Typical R-parity violating decay chain at the LHC. The primary vertex and the secondary 
vertices are highlighted. 



one of the main features searched for in usual searches for new physics. Furthermore the 
channels labeled as opposite sign could be found in usual supersymmetry (SUSY) searches 
as they include a considerable amount of ^7^, many jets and one isolated lepton pair with 
different signs. However, some searches remove events with muon pairs having invariant 
mass around the Z pole in order to dispose of Drell-Yan Z/7* — )■ // processes. Note that in 
the model presented in this work this cut would lead to a suppression of the signal. Other 
leptonic and semi-leptonic channels also contain opposite-sign lepton pairs but only small 
amount of $rp and therefore they are not considered in the usual searches, (c.f. [22,23]). 
Neutralino decays lead also to signatures containing same-sign lepton pairs but since no 
$'j' is present in these channels they are usually discarded in order to suppress various 
backgrounds [24]. 

If the value of C, is rather small one of the neutralinos will decay outside of the detector 
leading to signatures with large amount of $rp as shown in Table 2. The leptonic decays 
of one of the neutralinos inside the detector lead to a perfect opposite-sign signature. As 
mentioned above this signature can be hidden if one rejects events where the invariant 
mass distribution of the lepton pair is in the range of the Z-boson mass. Another strategy 
is the search for single lepton events with large amount of missing transverse energy. 

Thus the applicability and the reach of the usual SUSY searches applied to the model 
presented in this work depends crucially on the size of R-parity breaking. In order to 
further evaluate this statement we investigated a number of characteristic variables in 
supersymmetric events. The events were generated with PYTHIA as described in the next 
section, with the mSUGRA boundary conditions mi/2 = ttiq = 270, tan/3 = 10, a^ = 0, 
and n > 0. R-parity violating neutralino decays were taken into account. 

Fig. 5 shows the distribution of the P'y factors of the neutralinos. This factor enters 
the formula for the neutralino decay length and one sees from the plot that analytic results 
in the literature, which have been computed with P'y = 1, are correct within one order of 



category 



Xi decays 



LHC signature 



leptonic 



W+W-U -^ llllvv 

w+w+a -^ Tiiiuu 

W-W'U ^ llUvv 

ZW'lu -^ llUuu 

ZW+lu -^ lUluu 

ZZvv — > llllvv 



3j + 21 + 21 + $^ 



(opposite sign, 
$rp from, Z) 



ZW+lv -^ vvllyv 

ZW^lv — > vvllvv 

ZZvv — ;> vvllvv 



3j + ll + ll + :^T 



semi-leptonic 



W+W-ll -> jjlllv 

W+W+ll -^ jjlllv 

ZW+lv -^ lljjlv 



bj + 2l + ll + $T 



W+W~ll_ -^ j filly 

W-W'U ^ jjlUy 

ZW-lv -^ lljjlv 



5j + ll + 2l + $T 





ZW+lv -^ jjllvv 

ZW'lv -^ jjllvv 

ZZvv — ;> jjllw 


5j + ll + ll + $T 




W+W-ll ^ jjjjll 


7j + 11 + 11 


(same sign, 
no $t) 


W+W+ll -^ jjjjll 

w-w-U^jjjjU 


7 J + 21 
7 J + 21 


single lepton 
($T from Z) 


ZW+ly -^ jjjjly 
ZW'lv -^ jjjjlv 


7 J + ll + $T 
7j + ll + $T 




ZW^lv — )• vvjjlv 
ZW^lv —^ vvjjlv 


5j + ll + $T 
5j + ll + $T 


all-hadronic 
($T from Z) 


ZZvv -^ jjjjvv 
ZZvv — > vvjjvv 


7 J + $T 
5j + $T 


invisible 

($T from 2 Z) 


ZZvv — > www 


3j + $T 



Table 1: Possible final states if both neutralinos decay inside the tracking volume. 




Figure 5: P^ distribution of neutralinos at generator level for benchmark point HH27 (see Table 4)- 
The number of neutralinos corresponds to twice the number of the events scaled to 10fb~^ at 
^ = 7TeV. 




Figure 6: p distribution at generator level for benchmark point HH27 (see Table 4) and different 
values of the R-parity breaking parameter (. Generator level |zl„ is defined as sum over the px 
of i) neutralinos decaying outside of the detector (see Section 4-4) o,nd ii) all neutrinos produced 
inside of the detector. The number of events is scaled to 10 fb""*^ at ^/s = 7TeV. 
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category 


Xi decays 


LHC signature 


leptonic 

(opposite sign) 


W'l -> llu 
Zv -^ llv 


3j + ll + ll + $T 


single lepton 


wn -^ jji 
w-i^jfi 


5j + ll + $T 
5j + ll + $T 


all-hadronic 


Zv -^ jjv 


5j + $T 


invisible 


Zv — ^ vvv 


Sj + $T 



Table 2: Possible final states if one of the neutralinos decays outside the tracking volume. 

magnitude. The most important kinematic property connected with the neutralino decay 
length is the amount of missing transverse momentum p which is shown in Fig. 6 for 
different values of the R-parity violation parameter (. The missing transverse momentum 
was computed as the sum of the transverse momenta of all neutrinos produced in the 
detector before the hadronic calorimeter (r < 1800 mm, \z\ < 3700mm) and the transverse 
momenta of the neutralinos decaying outside the hadronic calorimeter. The p distribution 
of the R-parity conserving model ( = cannot be distinguished from the model with 
( = 1 X 10^^. However, the distribution is significantly different for ( = 3 x 10~^ since 
in this case most events have only very little missing transverse momentum due to early 
neutralino decays. This suggests that our model could only hardly be discovered in usual 
searches relying on ^x- A further analysis with full detector simulation is needed in order 
to properly evaluate the discovery potential of usual SUSY searches. 

Another general feature of models with relatively large ( is the large possible number 
of leptons in the final state, illustrated in Fig. 7. The generator level particles selected for 
this plot had to fulfill the criteria shown in Table 3 imposed in order to select leptons from 
hard processes which could be reconstructed in a realistic detector. The cuts on the vertex 
position represent a pessimistic estimate of the reconstruction efficiency (see Section 4.5). 

3.2 Search strategies 

As mentioned in the previous section one of the striking features of the presented model are 
events with secondary vertices and possibly many leptons in the final state. The search for 
a secondary vertex is crucial in order to ensure the R-parity violating nature of the decays. 
Possible search strategies can be optimized in order to find some of the channels described in 
Tables 1 and 2. It is remarkable that many channels allow for the full reconstruction of the 
neutralino mass: all decay chains including Z-bosons or hadronically decaying W^-bosons. 
The reconstruction of the neutralino mass from the particles produced in the Z-boson decay 
depends crucially on the full reconstruction of the secondary vertex, which is beyond the 
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Figure 7: The number of generated particles per event after selection cuts described in Table 3. 
The color code for the curves in both plots is given in Fig. 7b. The number of events is scaled to 
lOfb-i at ^ = 7TeV. 



scope of this work'^. This method of neutrahno mass reconstruction works also in R-parity 
conserving models where the neutralino decays into Z-boson and gravitino [17]. 

For example, one promising search strategy working for all ( values considered in this 
work is based on single lepton events with some number of hard jets and missing transverse 
energy larger than 90 GeV. After the preselection one could look for events where the lepton 
is coming from a secondary vertex and try to reconstruct the VT-boson mass from a jet 
pair. In the final step one could try to reconstruct the neutralino mass from the jets 
selected in the previous step and the lepton. However such study depends crucially on the 
knowledge of the detector response in the case of late decaying particles. A neutralino can 
decay in various detector components and lead to unusual signals. Furthermore the mass 
resolution is limited by the uncertainty in the jet energy scale and by the uncertainty in 
the determination of the jet momentum direction. 

We will focus our study on leptonic final states, which have a particularly clean sig- 
nature, and reconstruct the Z-boson coming from a secondary vertex. We will use only 
muon and track objects for which we assume to have modeled a realistic detector response 
(see Section 4.5). A possible background for this search are cosmic muons leaving no track 
in the detector. It is important to note that one would miss the signal in this channel 
entirely if one imposes a cosmic muon veto which rejects all events with muon pairs having 
no associated tracks (c.f. [22]). 



■^The four-vector pointing to the secondary vertex and the three- moment a of the leptons or jets from 
the Z-boson provide sufficient information for the reconstruction of the neutrahno mass. 
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particle transverse momentum pseudorapidity vertex position 

electron px > TGeV \ri\ < 2.5 r < 400 mm \z\ < 1300 mm 

muon pt > 6 GeV |?7| < 2.5 r < 4000 mm \z\ < 6000 mm 

Table 3: Cuts for the generator level particle selection for the study of particle multiplicity. 

4 Simulation of signal and background 

In this section we define a set of representative points in the parameter space of our model 
and describe the generation of the signal and dominant standard model (SM) background 
samples. In particular we examine the simulation of detector effects using the generic 
detector simulation DELPHES 1.8 [20] on signal and background in the presence of secondary 
vertices. 

4.1 Benchmark points 

A typical set of boundary conditions for the supersymmetry breaking parameters of the 
MSSM at the grand unification scale is given by equal scalar and gaugino masses, rriQ = 
mi/2. These boundary conditions lead to a bino-like neutralino Xi ^■s NLSP. We choose a 
representative value of tan/3 and set the scalar trilinear couplings to zero, 

ao = 0, tan/3 = 10. (24) 

Thus the universal gaugino mass remains the only independent supersymmetry breaking 
parameter which will be varied in the present study. 

Electroweak precision tests lead to lower bounds on the supersymmetric particle spec- 
trum. The LEP lower bound for the lightest Higgs at tan/3 = 10 is mn ^ 95 GeV 
[25]. In the present study the lightest superparticle spectrum corresponds to the choice 
rriQ = mi/2 = 270 GeV (HII27). At this benchmark point the NLSP is a neutralino with 
mass m^o = 105.8 GeV and the lightest Higgs boson has a mass ruh = 110.4 GeV. 

In order to probe the region of gluino and squark masses accessible at the LHC [26] we 
increase the gaugino mass parameter in four steps: mi/a = 350, 500, 650, 800 GeV. Some 
particle masses at these points are shown in Table 4. For the different benchmark points 
the production cross sections, calculated with PROSPINO 2.1 [27] at y/s = 7TeV, are given 
in Table 5. 

For the R-parity breaking parameter ( we choose the following values: ( = 3x 10^^ [10], 
C = 2 X 10-^ [14], 1 X 10-^ 5 X 10-9, 1 X 10-9, 5 x 10"^° and 1 x 10"^°. Note that for 
gravitino masses ms/a = (9(100 GeV) one has C ^ 1 ^ 10^^ [10]- 

4.2 Major backgrounds 

Neutralino decays always have W- and Z-bosons in the final state (c.f. Fig. 2). In our 
study we focus on the reconstruction of Z-boson decays to muon pairs. Therefore we only 
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GUT 


masses 




particle 


masses 






mo 


nil/ 2 


^xj 


rrih 


m~ 


"^« 


HH27 


270 


270 


105.8 


110.5 


662.4 


653.4 


HH35 


350 


350 


140.5 


112.5 


841.7 


831.8 


HH50 


500 


500 


205.7 


115.1 


1170 


1160 


HH65 


650 


650 


271.5 


116.7 


1492 


1481 


HH80 


800 


800 


337.8 


117.9 


1809 


1798 



Table 4: Definition of the benchmark points together with some particle masses; all masses are 
in GeV. 

consider SM backgrounds which lead to at least two muons in the final state originating 
from W- or Z-bosons: 

• tt production: W^-bosons from top quark decays. 

• Z production 

• Di-boson production {WW, WZ, ZZ) 

. Tri-boson production {WWW, WW Z , ZZW, ZZZ) 

Table 6a gives an overview of the background samples used in our analysis. We assume that 
pure QCD background can be efficiently suppressed in multi-lepton final states with high 
transverse momentum, particularly after imposing lepton isolation criteria (c.f. [23,28]). 

4.3 Event simulation 

All Monte Carlo samples were generated using parton distribution functions given by 
CTEQ6L1 [33]. For the simulation of the background we used MADGRAPH 4.4.44 [34] inter- 
faced with PYTHIA 6.4.22 [35]. 

Our simulation of the signal events relied on the following procedure. First, super- 
symmetric mass spectra were calculated with a modified version of SOFTSUSY 3.1.5 [36] 
assuming mSUGRA boundary conditions and R-parity conservation. The latter assump- 
tion is justified due to the tiny amount of R-parity breaking in our model. The SOFTSUSY 
version was modified in order to produce additionally to the spectrum the R-parity vi- 
olating neutralino decay width and branching ratios according to Equation (17). The 
SOFTSUSY mass spectra were fed into SDECAY [37] via the MADGRAPH homepage [38] in order 
to calculate the decay widths of the SUSY particles (besides the neutralino LSP). In the 
next step neutralino decay information was included into the SDECAY output. The signal 
process (production of gg, gq, qq and qq) was simulated with MADGRAPH and then given to 
PYTHIA for computation of all subsequent decays according to the SDECAY output as well 
as for parton showering and hadronization. Table 6b shows the signal samples used in our 
analysis. 
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(r{qq) 


a{qq) 


(^{gg) 


CT(tot) [fb; 
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1090 


682 


256 


208 


2236 


(739) 


(570) 
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(1566) 
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149 


38 
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385 


(105) 
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(25.2) 
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(265) 
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11.8 
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23.36 
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(10.1) 


(1.02) 


(0.206) 


(15.7) 


HH65 


0.579 


1.01 


0.0943 


0.0466 


1.73 


(0.216) 


(0.877) 


(0.0458) 


(6.37 X 10-3) 


(1.145) 


HH80 


0.0379 


0.0805 


5.37 X 10-3 


2.44 X 10-3 


0.126 


(0.0109) 


(0.0723) 


(1.98 X 10-3) 


(0.203 X 10-3) 


(0.0854) 



Table 5: Production cross sections at NLO (LO) at the benchmark points calculated with 
PROSPINO. 



The generic detector simulation DELPHES, tuned to the CMS detector, was used in order 
to account for effects of event reconstruction at the detector level. However, DELPHES 
describes the detector geometry solely in terms of angular variables, i.e. the detector 
is stretched infinitely in the radial direction. This approximation is sufficient for most 
studies involving prompt decays but is untenable in the case of late decaying particles. 
We overcome this obstacle by adding vertex information from particles at the generator 
level to objects at the detector level. Usually, this information is provided by the detector 
simulation. Our procedure is described in detail in the following section. We emphasize 
that a full detector simulation, which includes vertex reconstruction, needs to be done to 
improve our analysis. 

4.4 Muon reconstruction process 

Particles produced in the late decay of the neutralino will not be properly reconstructed in 
a real detector if the position of their vertex is beyond or even within the crucial detector 
component responsible for the respective identification. For example, an electron produced 
inside of the electromagnetic calorimeter will leave no track in the tracker and will therefore 
be identified as a photon or jet. In order to simulate the detector response to such events 
we use a detector geometry in the (r, z) coordinates, which is inspired by the CMS detector 
at the LHC (see Fig. 8). The angular position of the detector components is given by the 
CMS tune of DELPHES. 

In order to be as conservative as possible we only use muon and track objects for the 
present analysis, since these objects allow a simple simulation of detection efficiency losses 
due to the finite size of the detector. Namely, we assume that a muon can be reconstructed 
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Sample 


a [pb 


events 


ti 


163 


200000 


Z 


977 


700000 


w+w- 


47 


50000 


zz 


6.46 


50000 


w+z 


11.88 


50000 


w-z 


6.69 


50000 


w+w-z 


0.182 


15000 


w+zz 


0.040 


15000 


W+W-W+ 


0.146 


15000 


zzz 


0.015 


5629 



events 



HH27 


^ 1 X 10-9 
1 X 10-9 


22280 
222800 


HH35 


^ 1 X 10-9 

1 X 10-9 


10000 
100000 


HH50 


^ 1 X 10-9 
1 X 10-9 


10000 
100000 


HH65 


^ 1 X 10-9 
1 X 10-9 


10000 
100000 


HH80 


allC 


10000 



(a) Samples of SM background. The 
production cross sections are taken 
from [29-32]. 



(b) Samples of signal events for dif- 
ferent benchmark points (see Table 4) 
and ( ~ a X IQ-^ (a = 0.1, 0.5, 1,5, 

10, 20, 3o;. 



Table 6: Monte Carlo samples of SM background and signal events used for our analysis. 

as long as its vertex is in front of the muon chambers, and analogously a track can be 
reconstructed if it originates approximately in the first third of the tracker (This region 
is called pixel detector in Fig. 8). For the matching between generator level particles and 
objects reconstructed by DELPHES we use the distance in pseudorapidity f] and azimuthal 
angle 0, defined as AR = w (A0)2 + [Arjy. 

In the following we will call generator level muons, produced by PYTHIA, GenMuons, 
muons reconstructed initially by DELPHES muon candidates, and track objects reconstructed 
by DELPHES RecoTracks. Only GenMuons and Reco Tracks have the coordinates of their 
vertex. 

First, we perform the following p^ cuts on muon candidates and RecoTracks: 

. prifi) > 20 GeV, 

. pT(Track) > 15 GeV. 

These cuts are guided by our SUSY search strategy (c.f. Section 5), since we expect that 
muons coming from Z-boson decay have high px, and a sufficiently high pt cut can effec- 
tively suppress QCD fake leptons. Furthermore, DELPHES itself reconstructs only muons 
with Pt above 10 GeV. Additionally, these cuts were optimized in order to get a realistic 
muon reconstruction efficiency (see Section 4.5). 

In the second step vertex information is added to the muon candidates by matching 
with GenMuons: 

• A GenMuon is selected for matching with muon candidates if its vertex lies in front 
of the muon system : r^ = y/x'^ + y'^ < 4000 mm, |-2^| < 6000 mm (see Fig. 8). 
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Figure 8: Layout of one quarter of the generic detector used for particle identification. 

• The Ai? distance between each selected GenMuon and all muon candidates is com- 
puted. 

• A GenMuon vertex is added to the muon candidate closest in Ai?, if /S.R < 0.1 and 
GenMuon and muon candidate have the same charge. 

• Muon candidates with added vertex information are called RecoMuons. 

In the final step, muons with or without signal in the tracker are distinguished: 

• A Reco Track is selected for matching with RecoMuons if the track vertex lies in the 
following range: tt < 400mm , \zt\ < 1300mm. 

• Each selected RecoTrack is matched with the RecoMuon closest in AR, if AR < 0.1. 

• Matched RecoTracks and RecoMuons are called tracker muons. RecoMuons which 
cannot be matched with RecoTracks are called chamber muons. Each RecoMuon is 
therefore either a tracker muon or a chamber muon. 

After the reconstruction procedure one is left with two kinds of muon objects: (i) 
chamber muons which have no track in the tracker and are therefore reconstructed solely 
by the muon chambers, and (ii) tracker muons which have a track. The muon reconstruction 
process is depicted in Fig. 9. The AR matching condition has been optimized in order to 
get a realistic muon reconstruction efficiency (see next section). 
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Figure 9: Muon reconstruction process. 

4.5 Muon reconstruction efficiency 

In order to test our method of obtaining physically sensible objects we compute the muon 
reconstruction efficiency in the following way: 

• Muons are created as described above. 

• GenMuons are matched with RecoMuons without any constraints on the position of 
the GenMuon vertex. 

• The number of successfully matched objects is compared binwise (in bins of r and 
|2;|) with the number of all GenMuons. 

The second condition is necessary in order to see whether the assignment between Reco- 
Muons and GenMuons is correct. Since the matching procedure only relies on angular 
variables, it is possible that a RecoMuon originally matched with a GenMuon created in 
front of the muon chamber belongs in fact (i.e. has smaller angular distance) to a Gen- 
Muon coming from a decay inside the muon chamber or even outside of the detector. Such 
wrong matchings would be seen in the efficiency plot as efficiencies not equal to zero in 
regions where muons could not be detected by the detector defined above (r^ > 4000 mm, 
\z^\> 6000mm). 

Fig. 10 shows the computed muon efficiency in bins of r and \z\. As expected one 
sees a sharp decline in efficiency in the r plot at r^ = 4000 mm, where the hard cut 
applies. The decline in the z plot is gradually, since physical particles have to fulfill both 
r and \z\ criteria. The particles originating at small values of r and large values of |2;| are 
not reconstructed due to the limited pseudorapidity coverage of the muon detector. The 
efficiency stays at zero beyond r = 4000 mm and \z\ = 6000 mm as expected, confirming 
our method of muon reconstruction. We expect that the computed muon efficiency agrees 
within 15 % with efficiencies of present LHC detectors including losses due to muon-jet 
separation requirements. 
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Figure 10: Muon reconstruction efficiency for the benchmark point HH27. 

5 Search for the neutralino decay Xi ~^ ^^ 

As described in Section 3.2 our study is focused on the channel Xi ~^ Zu — >■ //"'"/i~i/. This 
channel possesses certain physical and technical advantages. On the physical side reliable 
muon identification is possible already in the early stage of the LHC data taking and one can 
assume that QCD background can hardly fake two muons at the same time. Furthermore 
this signal leads to spectacular events and has no easily identifiable SM background at 
all, as shown in this section. Additionally, the muon chamber is the detector component 
which is farthermost away from the primary vertex and hence one can expect that it will 
be possible to detect a significant number of clean late time decays even for very small 
R-parity breaking. On the technical side, muons seem to be the simplest objects for which 
a realistic detector response can be modeled within DELPHES (see Section 4.4), due to the 
limitations of this simulation in the presence of secondary vertices. 

The spectacular feature of this signal are opposite sign muon pairs with invariant mass 
close to the Z-boson mass, which have either associated tracks in the tracker with clearly 
visible secondary vertices or no associated tracks at all. Such muon pairs cannot be gener- 
ated by usual SM background as will be shown in the following. However, a similar signal 
can arise from cosmic muons traversing the detector. We could not create a Monte Carlo 
background sample for cosmic muons, and we simply assume that such background can 
be suppressed by use of the full timing information of the event: cosmic muons will first 
cause a signal in the muon chamber which is closest to the ceiling of the experimental hall 
followed by a signal in the opposite direction. 

An intrinsic background for the presented search are muon pairs from R-parity violating 
decays, where one muon is coming from the VT-boson decay while the other muon is coming 
either from the neutralino decay into the W^-boson in either of the two branches or from 
the W- or Z-boson decay in the second branch (c.f. Fig. 4). This background can be 
suppressed if one has access to the corresponding tracks by demanding that both of them 
originate from the same vertex. In the case of muons without tracks this background is 
irreducible. However it belongs itself to the signal one is looking for. 
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5.1 Event selection 

In order to find the signal we now employ a series of simple cuts on the reconstructed 
objects (muons, tracker muons, and chamber muons). As we will see, already with an 
integrated luminosity of only 1 fb~ at ^/s = 7 TeV a discovery of the benchmark scenario 
HH27 with C = 3 x 10"^ is possible. 

First, we perform a selection cut on the number of muons in the event: 

• A^(muons) > 2. 

We define two event classes depending on the number of tracker muons: 

• Class 1: the event contains at least two tracker muons A^(tracker muons) > 2. 

• Class 2: otherwise. 

From the description of the signal presented above, we implement additionally two sets of 
cuts depending on the class of the event. The cuts for Class 1 events are: 

• All possible invariant masses of opposite sign tracker muons are computed. The 
event passes the cut if at least one invariant mass is in the range of the Z-boson 
mass: 80 GeV < M^+^- < 100 GeV. If the event contains more than one appropri- 
ate combination of the tracker muons then the muons from the combination with 
invariant mass closest to the Z-boson mass are selected for further analysis. 

• (i(Vertex) > 5 mm: Each of the tracks associated with the two selected tracker muons 
should have a vertex which is further than 5 mm away from the primary vertex. This 
value is approximately one order of magnitude larger than the current resolution of 
the inner tracker (c.f. [28,39]). 

• A(i(Vertex)jj < 5 mm: The distance between the two track vertices should be less 
than 5 mm. 

• If the event fails one of the cuts it is classified as a Class 2 event. 
The cuts for Class 2 events are: 

• iV(chamber muons) > 2: If an event has less than two tracker muons than it should 
have at least two chamber muons. 

• All possible invariant masses of opposite sign chamber muons are computed. An 
event passes the cut if at least one invariant mass is in the range of the Z-boson 
mass: 80 GeV < M^+^- < 100 GeV. 

Since each Class 1 event is classified as a Class 2 event if it fails one of the cuts, no signal 
event is discarded because of the presence of muons with tracks not coming from neutralino 
decay. 
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2 
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Table 7: Cut flow for HH27 at ^/s = 7TeV. The number of signal events for Q 
(Q = 1 X 10^9 ) corresponds to an integrated luminosity o/lOfb-^ (^~ 100 £5-"*^^. 



3 X 10" 



Most events will fall into the second class. The analysis is then very simple and amounts 
to the search for events with muons without associated track in which the invariant mass 
of a muon pair lies in the Z-boson mass range. The cut flow is given in Table 7. As 
expected, no background events survived the cuts, since no standard model process should 
produce secondary vertices so far away from the primary interaction point. Although 
our background estimate has an uncertainty due to the limited statistics, we assume on 
physical grounds that no background events will pass the cuts if we increase the number 
of simulated events. However, the major uncertainty in this study, the number of the 
background events from cosmic muons, cannot be estimated with the present software. 
Therefore a full fledged analysis with full detector simulation which takes into account 
the cosmic muon background is needed. In the following we assume that this background 
can be efficiently suppressed with the full timing information of the event as described in 
the introduction to Section 5. Furthermore, we only estimate the systematic uncertainty 
due to the background and neglect statistical errors and the uncertainty of the muon 
reconstruction efficiency. 

The significance of the signal is computed with the profile likelihood method [40] in- 
corporated in the SIGCALC code [41]. We assume an integrated luminosity of lfb~^ at 
^/s = 7TeV LHC and a ten times higher Monte Carlo luminosity jCmc = Nh/(^b = 10fb~^ 
for all the background events. At this integrated luminosity 17 signal events and no back- 
ground events survive the cuts, which corresponds to a significance Zpi = 9.03. Instead, if 
one makes the pessimistic estimate that 1 background event from the cosmic muons passes 
the cuts one finds a significance Zpi = 6.39. Therefore we conclude that at the benchmark 
point HH27 with ( = 3x 10~^, R-parity breaking neutralino decays can be discovered with 
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Figure 11: Contour plot for the density of neutralino decays inside the detector per vaT^ ; the 
numbers on the horizontal boundaries of the detector components correspond to the total number 
of decays in the enclosed volume; mi/j = mo = 270 GeV, ^ = 3 x 10~^ and C = 10 fb""*^. 

the first inverse femtobarn of LHC data. 



5.2 Discovery reach at the LHC 

In tlie previous section we liave studied in detail tlie bencliniark point HH27: nii 



-/2 = mQ = 

270 GeV, wliicli yields tlie rather small superparticle masses m^o = 106 GeV, rUg ~ 660 GeV 
and rrLq ~ 650 GeV for the light quark flavors (cf. Table 4). From the decay rates given in 
Section 2.2 and the phase space factors shown in Fig. 3 one obtains for decay length and 
branching ratio into Z-boson/neutrino final states: 






4.5 m 



c 



10- 



BR{x° -^ Zu) ~ 0.14 



(25) 



Based on the production cross sections listed in Table 5 an integrated luminosity C = 
10 fb"^ yields about 22000 events and therefore 44000 NLSPs. 

We have studied this benchmark point for two different values of the R-parity breaking 
parameter: ^ = 3 x 10~^ and ^ = 1 x 10~^. For the larger value of ( one has cr^o ~ 50 cm. 
Hence, essentially all neutralinos decay inside the detector, most of them close to the origin. 
The spacial distribution of secondary vertices is displayed in the contour plot Fig. 11. Using 
BR{Z — )■ fi~^fi~) — 0.034 and the branching ratio given in Eq. (25), one concludes that 
there are about 200 events with a secondary x5"decay vertex, which contain a fi'^fi~ pair 



with M, 



/l^tM- 



Mz- This is consistent with the simulation which yields 174 events passing 
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Figure 12: Location of secondary vertices for the decays Xi ~^ -^^ with Z — t- /i^/^^ (blue crosses: 
inside pixel detector, black dots: outside pixel detector); the numbers on the horizontal boundaries 
of the detector components give the number of decays in the enclosed volume; rriiu = rriQ = 
270 GeV, C = 3 X 10"^ and C = 10 fb-^ 



all cuts (cf. Table 7). The locations of the secondary vertices of these events are shown in 
Fig. 12. 

For the smaller value of the R-parity breaking parameter, (^ = 1 x 10~^, the decay length 
increases to cr^o ~ 450 m. Now most neutralino NLSPs decay outside the detector. This 
is apparent from Fig. 13 where the total number of decays in the different subvolumina 
of the detector are given. Compared to C = 3 x 10~^, the number of decays inside the 
detector is smaller by a factor ~ 200, which roughly corresponds to the ratio of the decay 
lengths, as suggested in [15]. 

According to the simulation described in the previous section, for ^ = 1 x 10~^ an 
integrated luminosity of 100 fb~^ is needed to obtain 2 decays x? ~^ Zv — )■ /i+yU^z/, inside 
the detector. This number is very small compared to the total number of about 2000 
decays in the detector, which is a consequence of the tiny branching ratio into the chosen 
specific final state. It is likely that a substantially larger fraction of the events can be used 
in the search for a decaying neutralino. In [15] it has been argued that already 10 x\ decays 
inside the detector may be sufficient for the discovery of a decaying NLSP, which would 
require an integrated luminosity of only lfb~^. It remains to be seen whether for events 
with a secondary vertex and jets, signal and background can be sufficiently well separated. 

Let us now consider the benchmark point HH50: mya = wiq = 500 GeV, which implies 
the heavier superparticle masses m^o = 206 GeV and irtg ~ uiq ^ 1200 GeV for the light 
quark flavours (cf. Table 4). The phase space suppression is now negligible, fw — fz — ^, 
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Figure 13: Location of all neutralino decays inside of the detector (blue crosses: decays inside 
pixel detector; black dots: decays outside pixel detector); the numbers on the horizontal boundaries 
of the detector components correspond to the total number of decays in the enclosed volume; 
mi/2 = mo = 270 GeV, C = 1 x 10^^ and C = 10 fb-^. 
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Figure 14: Location of all neutralino decays inside of the detector (blue cross: decay inside pixel 
detector; black dots: decays outside pixel detector); the numbers on the horizontal boundaries 
of the detector components correspond to the total number of decays in the enclosed volume; 
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and one obtains for decay length and branching ratio into Z-boson/neutrino final states: 
cr^o~0.8m(^^j , BR{x1 ^ Zu) ^ 0.33 . (26) 

The total production cross section for these heavier gluino/squark pairs is about two orders 
of magnitude smaller (cf. Table 5), and therefore an integrated luminosity C = 10fb~^ only 
yields 460 NLSPs. 

We have studied this benchmark point again for the two different values of the R- 
parity breaking parameter ( = 3x 10~^ and ( = Ix 10^^. For the larger value of ( one has 
cr^o ^10 cm, and essentially all neutralinos decay inside the detector. The branching ratio 
into the considered final state is now somewhat larger, BR{xi — > Zu — > fi^fi'u) ~ 0.01, so 
that one expects about 4 events with this final state, which is consistent with our simulation. 
Hence, for this larger value of the R-parity breaking parameter and this benchmark point, 
the discovery of a decaying NLSP appears feasible already in the early phase of the LHC 

For ( = 1 X 10^^, the decay length is cr^o ~ 80 m and most neutralino NLSPs decay 
outside the detector. The spacial distribution of secondary vertices inside the detector, 
in total 12 for 10fb~^, is shown in Fig. 14. Due to the 1% branching ratio into the 
Zu — )• /i"^/i~z/ final state one then estimates that 1000 fb~^ will be needed for a discovery, 
which is consistent with our simulation. 

In Fig. 15 we have summarized the results of our simulations for the decay chain 
Xi — ^ Zu with Z — > fi'^fi". The benchmark points HH27-HH80 correspond to gluino 
and squark masses ranging from 650 GeV to 1800 GeV (cf. Table 4). The bands refiect 
the different number of events required for a 5cr discovery depending on the simulated 
background. The central value corresponds to 6 signal events (with luminosity C) with 
no background events for a simulated luminosity of 10 x C; the lower (upper) boundary 
represents 3 (13) signal events (with luminosity C) with no (1) background event for a 
simulated luminosity of 100 x C (10 x C). We conclude that with 10fb~^ a 5a discovery 
of a quasi-stable neutralino is possible for squark and gluino masses of 830 GeV (cf. HH35) 
and an R-parity breaking parameter (^ = 3 x 10~^, which is one order of magnitude smaller 
than the present astrophysical bound [10,14]. 

We expect that the sensitivity in the parameter ( can be significantly improved if 
also neutralino decays with jets are taken into account. Fig. 16 represents an estimate 
of the discovery reach for quasi-stable neutralino NLSPs at the LHC, assuming 10-20 
decays inside the detector (cf. [15]). The parameter space, which can be probed, is now 
significantly extended. As an example, with 10 fb^^ and squark and gluino masses of 
830 GeV (cf. HH35), one is now sensitive to C = 3x 10^^*^, which lies two orders of magnitude 
below the present astrophysical bound. Correspondingly, for heavier gluinos and squarks, 
rrig ~ rriq ~ 1480 GeV (cf. HH65), one can probe values of the R-parity breaking parameter 
down to C = 3 X 10"^. 
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Figure 15: 5a discovery reach in ^ for quasi-stable neutralino NLSPs via the decays x? — ^ Zu with 
Z — )• n~^fi~ . The different bench mark points correspond to gluino and squark masses between 
650 GeV and 1800 GeV; the bands represent different assumptions about the background (see text). 
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Figure 16: Estimate of the 5a discovery reach in C for quasi-stable neutralino NLSPs at the 
LHC; the lower (upper) boundary of the bands corresponds to 10 (20) decays inside the detector. 
The different bench mark points correspond to gluino and squark masses between 650 GeV and 
1800 GeV. 
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6 Summary and conclusion 

We have studied the decays of quasi-stable neutrahnos in supersymmetric extensions of 
the Standard Model with small R-parity and lepton number violating couplings which are 
consistent with primordial nucleosynthesis, thermal leptogenesis and gravitino dark mat- 
ter. As representative examples with neutralino NLSP we have considered five benchmark 
points, HH27-HH80, with gluino and squark masses ranging from 650 GeV to 1800 GeV. 

For the considered benchmark points the neutralino is bino-like, and the R-parity break- 
ing parameter (, which governs the neutralino decays Xi ~^ W^P and Xi ~^ Zu, is directly 
related to the partial gravitino decay width into photon and neutrino. 

The main goal of the present work has been to determine the range of the parameter 
( which can be probed at the LHC, for varying superparticle masses. As a conservative 
starting point, we have focused on events with a clean signature: cascade processes with 
jets where one of the produced neutralino NLSPs decays into Z-boson and neutrino, with 
a subsequent decay of the Z-boson into a muon pair. 

In Section 3 we have studied the qualitative signatures of these events, the /37 distribu- 
tion of the produced NLSPs, the jz* -spectrum and the number of leptons in the final state. 
A detailed simulation of signal and background events for the generic detector DELPHES 
(with CMS tune) has been described in Section 4, with emphasis on the reconstruction of 
muons. The crucial new element of our analysis has been the implementation of the finite 
NLSP decay length. 

The numerical results of our simulations have been summarized in Section 5. It is very 
instructive to look at the spacial distribution of decay vertices which have been given for 
different benchmark points and different values of (. One can clearly see that the sensitivity 
extends from decay lengths (9(50 cm), with NLSP decays mostly inside the detector, to 
values (9(500 m) where almost all NLSPs decay outside the detector. 

The results for the discovery reach for quasi-stable neutralino NLSPs roughly agree with 
the simple estimates which one obtains from the branching ratios into the Z{^^ ^~)v final 
state together with the assumption that these events are background free. It is remarkable 
that already with 10 fb~^ a 5cr discovery is possible for squark and gluino masses of 830 GeV 
and an R-parity breaking parameter (^ = 3 x 10~^, which is one order of magnitude smaller 
than the present astrophysical bound. 

It is likely that the severe restriction of our simulation to Z(yU+/i^)z/ final states can be 
relaxed and that a much larger fraction of events can be used for the analysis. Assuming 
optimistically that 10 decays inside the detector are sufficient for a discovery, one would 
be sensitive to values of C, down to 3 x 10^^° for squark and gluino masses of 830 GeV and 
10 fb~^. Alternatively, for gluino and squark masses of 1480 GeV one could probe values of 
the R-parity breaking parameter down to 3 x 10~^. For the same quark masses a luminosity 
of 100 fb~^ would improve the sensitivity in (^ by a factor of three. We conclude that for 
gluino and squark masses accessible at the LHC values of the R-parity breaking parameter 
can be probed which are far below the present upper bounds obtained from astrophysics 
and cosmology. 
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